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ABSTRACT 

We present Herschel SPIRE-FTS observations of Arp 220, a nearby ultra- 
luminous infrared galaxy. The FTS provides continuous spectral coverage from 



^Center for Astrophysics and Space Astronomy, University of Colorado, 1255 38th street, Boulder, CO 
80303 

^School of Physics & Astronomy, Cardiff University, Queens Buildings The Parade, Cardiff CF24 3AA, 
UK 

^ESA Astrophysics Missions Division, ESTEC, PO Box 299, 2200 AG Noordwijk, The Netherlands 
''Dept. of Physics & Astronomy, McMaster University, Hamilton, Ontario, L8S 4M1, Canada 
5JPL, Pasadena, CA 91109 

^Astrophysics Group, Imperial College, Blackett Laboratory, Prince Consort Road, London SW7 2AZ, 
UK 

^Blue Sky Spectroscopy Inc, Suite 9-740 4th Avenue South, Lethbridge, Alberta TIJ 0N9, Canada 

^UK ALMA Regional Centre Node, Jordell Bank Center for Astrophysics, School of Physics and Astron- 
omy, University of Manchester, Oxford Road, Manchester M13 9PL, U.K. 

^Istituto di Fisica dello Spazio Interplanetario, INAF, Via del Fosso del Cavaliere 100, L00133 Roma, 
Italy 

i°Dept. of Astronomy, University of Maryland, College Park, MD 20742 

"CEA, Laboratoire AIM, Irfu/SAp, Orme des Merisiers, F-91191 Gif-sur-Yvcttc, France 

^^Department of Physics & Astronomy, University of California, Irvine, CA 92697, USA 

^■^Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281 S9, B-9000 Gent, Belgium 

i^Laboratoire d'Astrophysique de Marseille, UMR6110 CNRS, 38 rue F. Joliot-Curie, F-13388 Marseille 
France 



- 2 - 



190 - 670 /xm, a wavelength region that is either very difficult to observe or 
completely inaccessible from the ground. The spectrum provides a good mea- 
surement of the continuum and detection of several molecular and atomic species. 
We detect luminous CO (J = 4-3 to 13-12) and water rotational transitions with 
comparable total luminosity ~ 2 x 10*^ L©; very high-J transitions of HCN (J = 
12-11 to 17-16) in absorption; strong absorption features of rare species such as 
0H+, H2O+, and HF; atomic lines of [C I] and [N II]. 

The modeling of the continuum shows that the dust is warm, with T = 66 K, 
and has an unusually large optical depth, with Tdust ~ 5 at 100 /im. The total 
far-infrared luminosity of Arp 220 is LpiR ~ 2 x 10^^ Lq. 

Non-LTE modeling of the extinction corrected CO rotational transitions 
shows that the spectral line energy distribution of CO is fit well by two tem- 
perature components: cold molecular gas at T ~ 50 K and warm molecular 
gas at T ~ 1350^^QQ K (the inferred temperatures are much lower if CO line 
fiuxes are not corrected for dust extinction). These two components arc not in 
pressure equilibrium. The mass of the warm gas is 10% of the cold gas, but 
it dominates the CO luminosity. The ratio of total CO luminosity to the total 
FIR luminosity is Lco/Lfir ~ 10~^ (the most luminous lines, such as J = 6-5, 
have Lco,j=6-5/-f'FiR ~ 10~^). The temperature of the warm gas is in excellent 
agreement with the observations of H2 rotational lines. At 1350 K, H2 dominates 
the cooling (~ 20 Lq/Mq) in the ISM compared to CO (~ 0.4 Lq/Mq). We 
have ruled out PDR, XDR and cosmic rays as likely sources of excitation of this 
warm molecular gas, and found that only a non-ionizing source can heat this gas; 
the mechanical energy from supernovae and stellar winds is able to satisfy the 
large energy budget of ~ 20 Lq/Mq. Analysis of the very high-J lines of HCN 
strongly indicates that they are solely populated by infrared pumping of photons 
at 14//m. This mechanism requires an intense radiation field with T > 350 K. 

We detect a massive molecular outflow in Arp 220 from the analysis of strong 
P Cygni line proflles observed in 0H+, H2O+, and H2O. The outflow has a mass 
> 10^ Mq and is bound to the nuclei with velocity < 250 km s"^. The large 
column densities observed for these molecular ions strongly favor the existence 
of an X-ray luminous AGN (lO"'"' ergs s^^) in Arp 220. 

Subject headings: galaxies: ISM, galaxies: starburst, ISM: molecules, techniques: 
spectroscopic, line: identification, molecular processes 



- 3 - 



Introduction 



Molecular gas is the raw material for the formation of stars within galaxies. However, the 
star formation process inputs energy and momentum back into this gas - and into the ISM in 
general - through radiation, stellar winds, and supernovae. The imp act of the feedb ack from 
both star formation and AGN activity on the ISM can be profound (jMaloneylll999l ). and the 
resulting change in the physical state of the gas must in turn affect the star formation process. 
In the general framework of galaxy evolution, disentangling the relative contributions of 
these processes for high-redshift objects, however, is an extremely difficult undertaking, 
due to their small angular size and faintness; in addition, many luminous high-^ sources 
are heavily obscured by dust, ruling out most of the traditional approaches for studies of 
galaxy energetics. Nearby galaxies such as M82, Arp 220 and NGC 1068, for which deep 
and relatively high resolution observations can be obtained at multiple wavelengths, are 
commonly used templates for calibrating the observations of high-z luminous dusty galaxies. 

The study of the star formation process and of the effect of AGN activity in nearby 
galaxies can be done in the far-infrared (FIR) and submillimeter (submm) spectral range, 
where a large variety of high-J rotational molecular transitions occur (CO, HON, HCO"'" etc.). 
These transitions have a large span in critical densities making them excellent tracers of the 
physical conditions of gas over a wide range in temperatures (10 - 1000 K) and densities (10^ 
- 10^ cm~'^). Also, their line ratios can be used as diagnostics for distinguishing between 
different energy sources responsible for the excitation of the gas, e.g., starbursts or AGN. 
The low-J transitions (up to J = 4-3) of these molecules trace the cooler gas and have been 
observed from ground-based radio and submm telescopes. These low-J transitions cannot 
distinguish between the AGN and starburst systems. However, the mid-J to higher-J tran- 
sitions, which can trace warmer an d denser gas and be used to test models that distinguish 



between AGN and starbursts (e.g. ISpaans fc MeijerinkI 120081 ) . are either diffi cult to observe 



or co mpletely inaccessible from the ground. Until the advent of Herschel (jPilbratt et al. 
2010l ) they had only been observed in a handful of extragalactic objects. 



The Very Nearby Galaxy Survey (VNGS) (PI: C. Wilson) is one of the guaranteed- 
time key projects of Herschel that is using both the Spectral and Photometric Imaging 



Receiver fSPI REl flGriffin et al. 



20101) . and the Photodetector Array Camera and Spectrom- 



eter (PACS) ( Poglitsch et al. 201ol ) instruments. This paper focuses only on the data from 



the SPIRE instrument, which is being used to obtain high S/N photometric and spectro- 
scopic maps of a sample of 13 nearby archetypical galaxies. The ima ging spectrometer on 
SPIRE is a Fourier Transform Spectrometer (FTS) ( iNaylor et al.ll2010l ). which provides con- 
tinuous spectral coverage from 190 - 670 fim (~ 450 - 1550 GHz). This region covers the 
spectral features of several key molecular and atomic species that are powerful diagnostic of 
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the physics and chemistry of the ISM. All the galaxies in this survey will have complemen- 
tary PACS observations either as a part of this program or another key project, allowing 
exploration of a full spectral view from 50 — 700/im. In this paper, we present the SPIRE 
observations of Arp 220. Until now only two nearby ga laxies observe d with SPIRE-FTS 



have been published. They are M82 ( jPanuzzo et al.ll2010l ) and Mrk 231 ( Ivan der Werf et al. 



20101). In both cases the mid- J to high- J CO ladder was detected. Modeling of this CO 
ladder showed that in M82 the excitation of molecular gas was governed by the starburst, 
whereas in Mrk 231 the excitation was dominated by the central luminous AGN. 

Arp 220 is the nearest Ultra Luminous Infrared galaxy (ULIRG) at a distance of about 
77 Mpc and z ~ 0.0181. It has Lpm ~ 10^^ L©, and is one of the most popular templates 
for studies of high-z du sty galaxies. It has two merging nuclei separated by ~ 1" (~ 360 pc 



at -DArp22o = 77 Mpc) (IScoville et al.l 119971 SC97 hereafter) and together they have a large 
reservoir of molecular gas (~ 10^° M©). Arp 220 has been observed extensively over the 
years across the electromagnetic spectrum. The extreme star formation environment in this 
galaxy provides an excellent laboratory to understand the processes affecting star formation 
and possibly AGN feedback. 

In this paper, we describe the state of the molecular gas by analyzing and modeling the 
new detections of high-J CO and HCN lines. We discuss sources of energy that could be 
responsible for the excitation of various molecular species detected in this spectrum. The 
issue of a hidden AGN has been long debated for Arp 220. We address this issue using 
the detection of rare molecular species that can only be detected by Herschel. We also 
present results from dust modeling by using the continuum measurements from SPIRE-FTS 
and SPIRE-photometer that adds data at wavelengths not accessible before. We summarize 
these results and discuss (in light of the new results) Arp 220 as template for broadband 
spectral energy distributions for high-z submm galaxies. 



2. Observations and Data Reduction 

Arp 220 was observed with the SPIRE-FTS in a high spectral resolution with FWHM ~ 
1.44 GHz (or ranging from 950 km to 280 km across the spectral range), single point- 
ing and sparse image sampling mode on 13 February, 2010. The total on-source integration 
time was 10,445 seconds. We used a deep dark sky observation with total integration time of 
13,320 seconds taken on 21 November, 2010 to account for the emission from the telescope. 
The FTS has two detectors arrays, called SLW and SSW, to cover the long and short wave- 
length ranges respectively, with a small overlap in wavelength. Arp 220 was also observed 
with the SPIRE photometer on 28 December, 2009 with a total on-source integration time 



- 5 - 



of 214s at 250 /iin, 350/im and 500 /iin. The calibration uncertainties in SPIRE-FTS are 
about 10% at the lower end of the frequency range but are as good as 5% across the rest 
of the spectrum. For the SPIRE photometer, the systematic uncertainty, which is related 
to the model for Neptune, is 5%. This is correlated among the three SPIRE bands. The 
random uncertainty, which is related to the ability of the instrument to repeat measurements 
of Neptune, is ~ 2%. This is independent for each of the SPIRE bands. 



We processed the SPIRE FTS data using the Herschel data processing pipeline ( iFulton et al 



20101 ). HIPE version 6.0. The spectrum is shown in Figured] and displays both emission and 
absorption features. Arp 220 is a point source for the diffraction limited SPIRE beam siz^, 
which ranges between FWHM of 17" — 40". The SLW and SSW continuum match perfectly 
across the overlap region. The flux densities from the SPIRE photometer were measured by 
fitting a two dimensional elliptical Gaussian function to the source in the timeline data. In 
this fit, the background was treated as a free parameter in the fit. The fit was performed to 
a inner region with radii of 22", 30", and 42" in the 250, 350, and 500 fim bands, respectively 
(which were selected because they contain the part of the beam profile that appears Gaussian 
and exclude the Airy rings) and using a background annulus with radii between 300" and 
350". The peak of the best fitting Gaussian function in each wave band corresponds to the 
flux density at that wavelength. The three SPIRE photometer measurements are overplotted 
in Figure [1] as red solid circles. They are in very good agreement with the FTS continuum 
within the statistical and calibration uncertainties of SPIRE-FTS and SPIRE-photometer. 

The spectral line shape of the FTS can be modeled to a good accuracy by a sine function 
with a FWHM of about 1.44 GHz. The ground-based observations measure large line velocity 



width s, of the order of 500 km s ^ in CO and HCN in Arp 220 (SC97) (also see lGreve et al. 



20091 ). Therefore, in our FTS spectrum we start resolving the broad spectral lines towards 
the higher frequency end. We use a sine function of variable FWHM to fit the emission 
and absorption lines using the MPFIT package in IDL. In Table 1 we report the integrated 
line fluxes in Jy km s~^ for emission lines and equivalent widths in fim for the absorption 
lines along with their 1— a statistical uncertainties. The observed spectral line shape suffers 
from small asymmetries from instrumental effects, and the line shapes could be affected by 
blending from weak lines. The effect on the line fluxes from asymmetries is expected to be 
< 2%. 



^The SPIRE beam shapes are not gaussian; the effective beam sohd angle can be found in the Herschel 
Observer's manual 
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Line Identifications 



The spectrum (Fig. 1) shows the FIR continuum and the detection of several key molec- 
ular and atomic species and Figure [2] shows zoom-in views of some of the spectral lines 
discussed in this section. In Table 1, we list their transitions, rest frequencies and integrated 
fluxes. The FWHMs mentioned in this section are from our sine line fitting. 

We detect a luminous CO emission ladder from J = 4-3 to J = 13-12 and several 
water transitions with total water luminosity comparable to CO. Comp ared to the velocit y 
widths measured in CO from the ground-based observations (e.g. SC97. iGreve et al.ll2009l ). 
the higher-J CO lines in our spectrum are much narrower. These high-J CO lines are not 
resolved and follows the resolution of FTS, which at the frequency of J = 12-11 line is ~ 310 
km s~^. This suggests that either the high-J CO emission is coming from only one of the 
nuclei or that it is a dynamically separate component from the low-J lines. On the other 
hand, the water lines are much broader with velocity FWHM of ~ 500 km s~^ implying that 
this emission could be coming from both nuclei. The observed velocity separation between 
the two nuclei in Arp 220 varies considerably from 50 km s~^ (Downes & Solomon 1998; 
Sakamoto et al 2008) to 250 km s~^(Scoville et al. 1997; Sakamoto et al 1999). Therefore, 
it is difficult to conclude whether the line widths are affected by the velocity shift between 
the two nuclei or by the dynamics of the molecular material. 

In addition to the emission line features, the spectrum also shows several absorption 
features. Most surprising are the detections of five very high-J HCN (J = 12-11 to J = 
17-16) lines in absorption; their low-J tr ansitions (J = 1-0 to J = 4-3) are detected in 
emission with ground-based observatories (iGreve et al.ll2009l . and references therein). Such 
high-J transitions of HCN have never been detected before in external galaxies, in emission 
or absorption. 

Strong absorption lines are detected from hydrides, including three transitions of OH"*" 
(see Figure [TU]) and one of CH+ (see Figure [2] (e)). We report four detections of 1120"*" 
transitions; two of them (1115, 1139 GHz) show st rong absorption and the other two, a t 
742 and 746 GHz (rest frequencies are reported in iBrudererl l2006l : lOssenkopf et al.ll2010l ). 
appear in emission. Before Herschel, H2O"'" was only detected in comets. The transitions at 
1115 and 1139 GHz (see Figur e [T0|) were first detecte d in Galactic mole cular clouds and in 



extra galactic objects by HIFI f lOssenkopf et al.ll2010l ) and SPIRE-FTS ( Ivan der Werf et al. 



20101). We have the first detection of transitions at 742 and 746 GHz (see Figure E^c)). 
The detections of OH"*" and H20^ are very important in Arp 220 for modeling the water 
transitions because they are the major intermediaries in the ion-neu tral chemistry network 
producing water in the ISM ( iGerin et al.ll2010l : iNeufeld et al.ll2010al ). and their detection is 
very difficult from current ground-based observatories. Herschel HIFI and SPIRE-FTS also 
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made the first detection of tlie HF 1232 GHz line (iNeufeld et al.ll2010bl : Ivan der Werf et al. 



20101 ). and we detect it in absorption in Arp 220. Tlie spectrum also sliows several features 
of Nitrogen hydrides including NH, NH2 and NH3. 

The rest frequencies of all the above detections of the molecular species are consistent 
with the systemic z = 0.0181 of Arp 220. However we detect two emission features at ~ 532.6 
and ~ 536.8 GHz (see Figure [2] (a)) that are close in frequency to the J = 6-5 transitions of 
HCN and HCO^; their low- J transitions (up to J = 4-3) have been seen in emission from the 
ground. These features also have larger FWHM (~ 1000 km s~^) compared to the spectral 
resolution at these frequencies. The redshifted frequencies of both HCN and HCO"'" are 
offset from the systemic velocity of Arp 220 by ~ —0.9 GHz (-500 km s~^) and ~ —1.7 GHz 
(-1000 km s~^), respectively. Such large shifts are not seen for any other detections in our 
spectrum. Since we see very high-J lines of HCN in absorption it is possible that the HCN 
J = 6-5 line is partially absorbed by the same gas producing HCN in absorption at higher 
frequencies. This effect could shift its apparent central frequency. The case for HCO"'" is 
even more problematic. The rest frequency of the emission feature is in perfect agreement 
with the rest frequency of HOC"*" molecule instead of HCO"*". Since the width of this feature 
is about ~ 1000 km s~^, these two molecules could be blended. However, identifying this 
feature as HOC"*" would imply that it has comparable to or higher line strength than HCO"*". 
This would be extremely unusual and has never been observed in other galaxies or in Galactic 
molecular clouds. All observations of HOC"*" and HCO"*" show that the latter is at least 100 



- 1000 times stronger (IMartin et al.ll2009l : ISavage fc Ziurysll2004l ). In Table 1 we mark these 



identifications with a question mark. 

Among the atomic species we detect [C I] and [N II] lines (see Figure [2] (d) & (i)). 
The detection of the two [C I] lines at 492 and 809 GHz have been rare in the past. The 
forbidden [N II] line at 1462 GHz has been inaccessible from ground and rarely accessible 
from space. The Herschel FTS has been very successful in detecting them in several Galactic 
and extragalactic targets. The [C I] (1-0) 492 GHz line survey for nearby galaxies was done 
by Gerin & Phillips (2000), in which this line was detected in Arp 220. We report the 
first detections of [C I] (2-1) and [N II] lines in Arp 220. The [N II] line width is broader 
(~ 400 km s~^) than expected from the FTS resolution of about 280 km s~^ at these high 
frequencies, implying that we might be resolving the [N II] emission from the two Arp 220 
nuclei. More discussion of these two atomic lines is presented in Section [3 
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4. Dust in Arp 220 

In addition to line fluxes, the FTS also provides a good measurement of the continuum 
at wavelengths that have not been accessible before. We model the dust spectral energy 
distribution (SED) by using the combination of FTS continuum data (listed in Table 2) 
and other data in the literature, providing wide wavelength coverage. The modeling yields 
constraints on dust temperature (T), spectral index and the frequency (uq) at which the 
dust optical depth becomes unity. 



We combined FTS continuum with ISO-LWS continuum data (IBrauher et al.ll2008l ) pro- 
viding a wavelength coverage from 50 — 670 /im. We chose ISO-LWS over other observations 
because it measures several data points covering the peak of the SED. In addition, the cal- 
ibration of the two data sets agree extremely well. We also added three continuum data 
points from the SPIRE photometer at 250 /im, 350 /im and 500 /xm. We fit these data with 
the following modified blackbody model: 



where 



& (T) 



hu 

QkT — \ 



(2) 



We generated a grid of models over a range in T, /3 and z/q. Using this grid we produced 
likelihood distributions of T, /3 and vq by comparing the model continuum fluxes to the 
observed ones (corrected for ZArp22o = 0.0181). In addition to statistical uncertainties we also 
included calibration uncertainties using a covariance matrix, a much better approach than 
adding them in quadrature. Figure [3] shows the highest likelihood model fit to the SPIRE- 
FTS (blue - open circles), SPIRE-photometer (red - open triangles; listed in Table 2) and ISO- 
LWS (red - open squares) data. The reduced chi-squared of this fit is 0.5, indicating that our 
calibration errors are likely too conservative. This excellent fit to the data is consistent with 
only a single temper ature component. We overplot data from IRAS in gree n-solid squares 
JSanders et al.lboosh . SC UBA in green - so hd triangles flPunne et allboooh . ISOPHOT in 
purple - open diamonds (IKlaas et al.l l200ll ) and Planck in orange open triangles (Planck 
Collaboration 2011, larXi v:1101.2041^. which also show good agreement with the fit down 
to 15 /im, below which the single temperature component fit breaks down. The contours in 
Figure H] show 68%, 95% and 99% confidence levels for the three model parameters. The 
addition of SPIRE data provides far better constraints on the dust parameters, as indicated 
by the compact contours, than was possible with previous measurement alone. The highest 
likelihood values and their l-o" ranges are: T = 66.7 K (66.4 - 66.9), f3 = 1.83 (1.82 - 1.85), 
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and uq = 1277 GHz (1262 - 1287). Ge nerally, the value of the emissivity index, /3, is fixed at 
2.0 for the sihcate and graphite dust (jPraine fc Leelll984j ). Treating /3 as a free parameter 
in our analysis gives a value of ~ 1.8, which is very close to the model value. The effect of 
dust extinction in Arp 220 is very large. The dust optical depth is unity at ~ 240 fim and 
~ 5.0 at 100 /im. The total hydrogen column density is ~ 10^^ cm~^ using the 250 /im dust 



cross-section from iLi fc Draind ( l200ll ). which has a systematic uncertainty close to a factor 
of 3. All the emission line fluxes in our spectrum need to be corrected for this large dust 
extinction. We estimate a dus t mass, Mnnst ~ 10 ^ Mq, for this model using the flux at 250 
/im and k,25o = 0.517 m^ kg~^ ( Li fc Draine 2001). This gives the gas-to-dust mass ratic[§ of 



100 using the molecular gas mass of 
fll997h . 



10^° Mq for Arp 220 measured by iScoville et al. 



The large dust optical depths at submm and FIR wavelengths in Arp 220 are not unex- 
pected. There is a high column density of molecular gas measured from CO line observations 
(e.g. SC97). The FIR spectrurn shows both [01] 63 /im and [N II] 122 /im lines in absorption 
( iGonzalez- Alfonso et al.l 120041 : iFischerl 120001 ) , compared to other galaxies where these lines 
usually appear in emission, suggesting a stronger dust continuum in Arp 220. The [C II] 
158 fim line deficit relative to the FIR continuum (Lc,+ / Lp -m ~ 2 x 10~'^) in th is galaxy is 



one of the most extreme compared to other ULIRGs (e.g.. iLuhman et al.ll2003[ ). From an 
interferometric map of Arp 220 at 860 /im, Sakamoto et al] ~( 2008 ) estimated a large dust 
optical depth, of ~ 1, even at this wavelength. In addition, if the observed CO ladder is not 
corrected for the large dust extinction (see section , then the temperature derived for the 
molecular gas traced by CO from our radiative transfer modeling will be in large disagree- 
ment with the temperature derived from the mid-IR H2 lines observed from Spitzer. Several 



previous works have argued for large dust optical depths in this galaxy fe.g.. iLisenfeld et aL 



2000; 



Fischer 



2000 



Papadopoulos et al. 



K!laas et al. 2001; Gonzalez-Alfonso et al. 2004; Downes fc Eckart 



2007 



2010h . Our result of large dust optical depth is consistent with all these 



studies. However, the addition of continuum data from SPIRE has reduced the uncertainties 
on the dust parameters by an order of magnitude compared to the previous studies. 

Multiple studies have shown that dust SEDs for ULIRGs may not only be fit by a 
single temperature optically-thick dust component but also by multiple optically-thin (as- 



(Klaas et al. 


2001; 


Blain et al. 


2003) 



temperature component model to the Arp 220 continuum data did not produce realistic 



^The molecular to atomic gas mass fraction is much higher for ULIRGs (IMirabel fc Sanderslll989l ) com- 
pared to quiescent galaxies, so we ignore the contribution of the atomic gas while calculating the gas-to-dust 
mass ratio. 
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results. Our optically-thin 2-component model fit gives ~ 25 K and ~ 46 K for the cold 
and warm temperature components, and beta = 2. Even though these temperatures are 
realistic the dust mass for this mo del is not: i t yiel ds an estimated dust mass of ~ 10^ 



Mq, using the flux at 250 fim. A iDraine fc Lil (120071 ) model fit for an optically thin dust 
assumption also produces a large dust mass (~ 6 x 10^ M©). This is an order of magnitude 
larger than the optically thick result of ~ 10^ Mq. For M(H2) ~ 10^° in Arp 220 (SC97), 
the gas-to-dust mass ratio is only 10 in the optically thin case. Using our preferred value of 
M(H2) = 5.2 X 10^ Mq derived from modeling of the CO line energy distribution (see section 
5 and table 3), the gas-to-dust mass ratio is still only ~ 5 in the optically thin case; using 
larger value estimated by SC97 would only raise it by a factor of 2. Even after considering 
a factor of three uncertainty in K.2r,n, this ra t io is extremely low (expected to be > 100, e.g., 
Hildebrandl Il983l : iDevereux fc Young! Il990l : ISanders et al.lll99lh . Therefore, modeling the 
dust with a general modified blackbody model in which no prior assumption is made about 
dust optical depth is strongly favored. 

Historically, fitting a single thermal component to model dust emission between 24 
lim and 600 fim has been considered unrealistic. This is because the dust is expected to 
be distributed at different distances from the heating source; the dust observed at 60 /xm 
would be much closer to the heating source than the dust observed at 500 fim. Also, even 
for dust heated by a uniform radiation field, transiently-heated smaller dust grains will 
typically produce emission that peaks at shorter wavelengt hs, while the larger dust grains 



will a ppear to be in thermal equilibrium at 15-30 K (e.g., iLi &: Drainell200ll : iDraine &: Li 
20071 ). However, in case of Arp 220, the dust is optically thick at these shorter wavelengths, 
and so transiently-heated small grains as well as grains closer to the heating source are not 
visible. From the average dust measurements we find that the FIR luminosity is dominated 
by a skin at ~ 66 K. 



We estimate a total FIR luminosity LpiR, 



1.77 X 10^^ Lq by integrating the 



■L(S-lOOOnm) 

model SED. Using this value of LpiR and T = 66 K in the Stefan-Boltzmann relation we 
obtain a source size of about 240 pc radius (~ 0.64" using -DArp22o = 77 Mpc). This gives 
a source solid angle of 4.1 x 10~^^ sr, which is a factor of ~ 1.6 smaller than measured for 
the CO (1-0) hne by SC97. Our single temperature fit suggests that the dust emission is 
coming from one nucleus alone. This is also supported by our high-J CO line observations, 
which have much smaller line velocity widths compared to the low-J CO lines, implying 
that the high-J CO emission is predominantly coming from one nucleus. This conclus ion is 
inconsiste nt with the submm inte rferometric observations from (jSakamoto et al.ll2008l ) (860 
/im) and (IMatsushita et al.l l2009l ) (435 /im). These studies suggest that the dust emission 
in the western nucleus is more compact than the eastern nucleus and the luminosity of the 
western nucleus is either comparable to or higher (2-3 times) than the eastern nucleus. The 
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sizes and the emission of the two nuclei are hkely to change with wavelength such that this 
trend may not hold over the wavelength range that we are observing. 



5. CO: Spectral Line Energy Distribution and Non-LTE modeling 



The CO line luminosities are affected by the large dust optical depths derived in the 
previous section. We correct for this dust extinction by assuming that the dust and gas are 
well mixed; i.e., / = /o(l — e~^^)/rA (IHaas et al.ll200ll . and references therein). Figure |5] (left 
panel) shows the observed extinction corrected luminosity distribution of the CO ladder from 
J = 1-0 to J = 13-12. The extinction correction factor {Iq/I) ranges from 1.08 at 450 GHz to 
1.95 at 1600 GHz. The solid circles are the FTS measurements from this work and the blue 



triangles are the average of the ground-based observation s compiled bvlCreve et al.l ( l2009l ). 
In addition, we show new observations from the JCMT (jPapadopoulos et al.ll2010l ) of CO 
(6-5) line and from Z-Spec (private comm.) of CO (2-1) line for comparison. The JCMT flux 
is lower by a factor of 1.3, likely missing some flux owing to large velocity widths in Arp 220. 
The Z-s pec CO (2-1) meas urement is higher than all the other ground-based measurements 
listed in lGreve et al.l (120091) but in agreena e nt wi th the JCMT CO (2-1) measurement within 



their 1—a uncertainties. 



Matsushita et al.l (120091 ) obtained an interferometric map of CO J = 
6-5 line using SMA. Figure 7 of their paper indicates that the CO J = 6-5 emission is mainly 
coming from only one of the nuclei compared to the low-J CO transitions (J = 1-0 to J = 
3-2), which clearly show emission from both nuclei. This is consiste nt with our observations 
of the narrower line widths of the high- J CO transitions. However, [Matsushita et al.l (120091 ) 
are missing ~ 70% of the flux in this line, obtaining integrated line flux of only 1250 ± 250 
Jy km s~^. 

The total CO luminosity is about 2 x 10*^ Lq and is dominated by the mid-J to high-J 
CO transitions. The shape of the CO spectral line energy distribution (SLED) in Arp 220 is 
similar to that of M82 (PIO), a starburst galaxy, in which the CO line fluxes are highest for 
the mid-J lines and then fall off at higher-J. This is in contrast to Mrk 231 (Van der Werf 
et al. 2010), in which the CO SLED rises up to J = 5-4 and remains flat for the higher-J 
transitions, which can be explained by the presence of a central X-ray source illuminating 
the circumnuclear region. The shape of the CO SLED has been sug gested as a discriminant 
betw een star formation dominated versus XDR dominated systems (Ivan der Werf fc Spaans 
20091 ). It will be possible to test this in the near future by compiling CO SLEDs for a sample 
of nearby galaxies observed by Herschel from several key projects. However, a practical 
application of using the shape of the CO SLED as a discriminant may be complicated in 
systems where the luminosities of star formation and AGN are comparable. In Arp 220, the 
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shape of the CO SLED is consistent with star formation but we show later (in Section 8) 
that there is also evidence for a luminous AGN from observations of other molecules in the 
FTS spectrum. This is discussed further in Section 9. 



5.1. Non-LTE Modeling 

We used the full observed CO SLED from J = 1-0 to J = 13-12, with the exception of the 
10-9 line, which is blended with a water line (see Figure |2] (f)), to mo del the CO excita- 



tion a nd radiative transfer using a variant of the non-LTE code, RADEX (Ivan der Tak et al. 



20071 ). This code computes the intensities of molecular lines by iteratively solving for sta- 
tistical equilibrium, using an escape probability formalism. In our case we used the escape 
probability formalism for an expanding spherical shell (large velocity gradient, LVclfl). The 
code starts from an optically thin case to generate the initial guess for the level populations, 
then iterates until a self consistent solution is achieved such that the optical depths of the 
lines are stable from one iteration to the next. The inputs to RADEX are the gas den- 
sity (^Ha), the kinetic temperature (Tkin), and the CO column density per unit line width 
(Nco/dv), which sets the optical depth scale. 

We used this code to compute CO intensities for a large grid in Tkin (10-3000 K), n^^ 
(10^-10^ cm~^), and Nco/dv (lO^^'^-lO^^-^ cm~^ (km s~^)~^). From this grid, we gener- 



ated likelihood distributions, adapting the method described in IWard et al.l ( 120031 ) (also see 
Kamenetzky et al. 2011), for Tkin, n]i2, and A^^co, by comparing the RADEX and observed 
line fluxes. 

To avoid any non-physical solutions we applied three priors in the likelihood analysis. 
The first one limits the ^^CO column density such that the total mass of the molecular gas 
producing the CO lines cannot exceed the dynamical mass of the galaxy according to the 
following relation: 

— — < — — = 1.6 X 10 cm (kms ) (3) 

dv fJ^rnii^ AAV 

where the dynamical mass of the disk M^yn = 1 x 10^° Mq (SC97), the abundance of CO 
relative to H2, xqo = 3 x 10"^, A is the area in cm^ covered by the source and is given by A 
= (Angular diameter distance)^ xfi^, and fi = 1.4 is the mean molecular weight in units of 
mn- The source size, Qg, is different for the high-J and the low-J lines; the values we used 



■^The form of escape probability for LVG and uniform sphere models are very similar and hence either 
choice will not introduce any significant difference in the CO flux values. 
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are described later in this section. To get the total CO column density, the Nco/dv should 
be multiplied by the line velocity width {AV). We use AV = 350 km s~^ for all the FTS 
lines, which comes from the line fitting of J = 12-11 (which lies at a frequency where the 
spectral lines are resolved). The second prior limits the column length to be less than that 
of the entire molecular region according to 



^^^<2i?d~ 400pc (4) 

where Rd ~ 200 pc is the radius of the emitting region. This prior will limit the density 
at the lower end and column density at the higher end. 

The third prior is the absolute surface brightness of the CO lines that is set by the choice 
of our source size and an area filling factor of unity. This prior will effect the column density 
such that it will systematically shift to higher values for smaller assumed source sizes. 

We first modeled the J = 4-3 to J = 13-12 transitions and found that the highest 
likelihood model, with a temperature of ~ 1350 K, provides a good fit to all the transitions 
with the exception of J = 4-3. The model fit to the observed line temperatures is shown in 
Figure [5] (right panel). The mid- J to high- J lines are tracing very warm molecular gas. The 
predicted Ray leigh- Jeans temperature of the J = 4-3 line is underestimated with respect 
to the observed temperatures by this warm component, implying a second temperature 
component. We then modeled the low-J lines from J = 1-0 to J = 4-3, after subtracting 
the contribution from the warm component, and found that these lines are fit well by a cold 
temperature component of ~ 50K as shown by the dashed line in the right panel of Figure 
[H The total of the two is shown by the dot-dashed line. 

The likelihood distributions of Ttin, nn^^ Nqq and thermal gas pressure (Tkin x ) '^^^ 
shown in Figure [6] for the warm (red) and cold (blue) components. The likelihoods for the 
warm component are narrow and provide very tight constraints on these parameters. Note 
that neither priors in Equation (3) and (4) were relevant for modeling the warm CO. In 
fact, only the third prior is important; the absolute line temperature set by the source size 
{Qg = 4.1 X 10~^^ sr) calculated in section HI This source size was calculated assuming an 
optically thick emission, which is a reasonable value to use in absence of spatially resolved 
CO maps at these high frequencies. As mentioned previously, any change in the source 
size will affect the likelihood distributions, particularly the column density, which will shift 
systematically to lower or higher values. The effect on density and temperature likelihood 
is minimal. The optical depths of mid-J to high-J transitions range from 3.3 to 0.1, peaking 
at J = 6-5. The likelihoods for the cold component are also very well constrained. For the 
modeling of the cold component we used the source size of 6.4 x 10^^^ sr and line velocity 
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width of 500 km as measured by SC97 for the CO J = 1-0 hne. In this case also the 
third prior is more important. The mass prior in Equation (3) does hmit the column density 
at the higher end but just barely. The low- J transitions are optically thick. 

The highest likehhood values of Tkin, ^Ha, -^co, pressure, mass, and luminosity are listed 
in Table 3 for the two components. Their pressure likelihoods do not overlap implying that 
these two components are not co-spatial. In addition, the observed velocity widths are also 
significantly different between these components. Even though the mass of the warm gas is 
only about 10% of the mass of the cold gas, the luminosity and hence the cooling is still 
dominated by the warm high-J CO lines. 

The molecular gas traced by CO in Arp 220 is similar to M82, where PIO also found 
warm and cold temperature components. However the temperature of these components were 
smaller in M82 by a factor of three and two for the warm and cold components, respectively. 
This is not unexpected given the higher FIR luminosity of Arp 220. 

The CO emission that we observe, in both the warm and cold components, arises in 
fairly lo w density (n ^ 10^ cm~^). Undoubtedly denser gas components are present; for 



example, iGreve et al.l ( 120091 ) argue for gas with n ~ 10^ cm~^ and T ~ 50 — 120 K based 
on their observations of CS lines. That we find no evidence for such a component indicates 
only that the contribution from this gas to the observed CO emission is small. 



5.2. Comparison with H2 Observations 



Higdon et al.l (120061 ) conducted a survey with the Spitzer Space Telescope of warm 



molecular gas in nearby ULIRGs (including Arp 220) by measuring the optically thin ro- 
tational lines of H2. They list line intensity measurements for the S(l), S(2), S(3) and 
S(7) lines, and provide upper limits on the S(0) line. We apphed a dust correctior0 to 
these line intensities using the mixed du st model arid the mid infra-red extinction curve of 
Kruegel fc Siebenmorgenl (119941 ) (also see lHaas et al.l (l200ll )). and used them t o calculate the 



excita tion temperature following the method described in sections 3 and 4 of iHigdon et al. 
(120061 ). We found the excitation temperature to be ~ 1380 K for the warm H2 gas using 



the S(3) and S(7) lines, in excellent agreement with our warm CO temperature. We note 
that without making the extinction correction for H2 and CO, the derived gas temperatures 
are grossly underestimated: ~ 260 K and ~ 400 K respectively. We estimated a total lu- 



"'We found the extinction correction factor, Iq/I ^ 3.7, for the S(l), S(2) and S(3) transitions, and ^ 5.3 
for the S(7) transition. 
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minosity from the observed H2 lines to be ~ 2 x 10^ Lq, and infer a warm H2 mass 
of ~ 3 X 10® Mq. Both the luminosity and mass of warm H2 are consistent with our 
warm CO measurements. This implies that the warm H2 gas observed by Spitzer is the same 
molecular gas traced by the high-J CO lines measured by Herschel. 

The total H2 cooling at 135 K is about 20 Lq/Mq using the theoretical H2 cooling 



curve of iLe Bourlot et al.l ( 119991 ) . This is much larger than the total cooling from CO of 
~ 0.4 L0/M0 implying that H2 is still the dominant coolant at high temperatures. The total 
H2 cooling in Arp 220 is almost an order of magnitude large than M82 (PIO); a consequence 
of the much higher temperature in Arp 220 than in M82 (1350 K vs 500 K). A cooling of 
20 Lq/Mq would imply that the warm molecular gas mass in Arp 220 would produce a 
luminosity of about 9 x 10^ Lq; a large number but still a small fraction (5 x 10"'^) of the 
total FIR luminosity of Arp 220. 



5.3. Heating Sources of the Warm CO Molecular Gas 

Here we consider four major sources of energy that could be responsible for producing 
the luminosity of the warm CO: heating by UV star-light in photo-dissociation regions, 
heating by X-rays in X-ray dominated regions, heating by cosmic rays (CR), and mechanical 
heating. 



5.3.1. Photo-Dissociation Regions 



Photo-dissociation regions are neutral regions present at the skins of molecular clouds 
that are illuminated by intense far-UV (FUV) radiation fields. This FUV radiation can 
strongly influence the gas chemistry in these regions and is also the most important source 
of heating. 



We explore a PDR model grid (IKaufman et al.l l2006l : IWolfire et al.l I2OIOI ) with predic- 
tions for CO line fluxes from J = 1-0 to J = 30-29, generated over a large range in density 
( 10^ - 10^ cm-3) and incident FUV flux. Log Gq (-0.5 - 6.5), where Gq is the FUV flux 
in units of the local interstellar value also sometimes referred to as th e Habing interstella r 
radiation field i.e., FUV flux = 1.3 x 10~^ x Go erg cm'^ s"^ sr"^ JWolfire et al.lboioh . 
Also, these models assume that FIR flux = 2x FUV flux. We first compare the ratio of the 
CO J = 6-5 line flux to the FIR surface brightness in Figure [7] (top left). The red dotted 
contour lines are the ratio of the model CO (6-5) line flux to FIR surface brightness, the 
green dashed line is the observed ratio, and the black dashed lines show the 1 — a density 
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limits from Table 3. The FIR flux is labeled on right hand axis of the left panel. The ob- 
served ratio can be produced within the acceptable density limits. However, the model FIR 
surface brightness is only ~ 0.1 erg cm~^ s~^ sr~^ compared to the observed FIR surface 
brightness, which is about ~ 70 erg cm~^ s~^ sr~^, 700 times larger. The PDR models fail to 
explain the observed FIR surface brightness. In principle, one can stack up multiple PDRs 
to reproduce the observed FIR. In our case, about 700 PDRs would have to be stacked to 
match the observed FIR surface brightness. However there is a limit to the number of PDRs 
that can be stacked, which is set by the characteristic column density of about 10^^ cm~^. 
The observed column density for warm molecular gas is (from Table 3) ~ 10^^ cm~^ (using 
^co = 3 X 10~^), allowing a maximum of 10 PDRs for stacking, not sufficient to match 
the observed FIR surface brightness. Moreover, in Figure [7] (top right) we compare the line 
ratio of CO J = 12-11 to CO J = 5-4. The observed line ratio (green line) of 0.7 cannot be 
produced by PDR models within the accepted density limits. In fact, the densities required 
to match the observed ratios would have to be greater than 10^'^ cm~^, completely excluded 
by our likehhood analysis. 

From these two comparisons, we conclude that PDRs cannot produce the observed warm 
CO luminosity and hence are not the dominant source of heating. They can only provide 
about 1.4% of the observed FIR flux. 



5.3.2. X-ray Dominated Regions 



The XDR modeling is done using an updated version of the code described in lMaloney et al 



(119961 . MHT). We assume a hard X-ray (1-100 keV) luminosity of 10^^ erg s~^, as would be 
expected for an object of Arp 220's luminosity if a large fraction of the bolometric luminosity 
(~ -Z^fir) were produced by an AGN. A power-law index a = 0.7 is assumed. Unlike the 
models in MHT, we assume a sharp lower energy cutoff of the incident spectrum at 1 keV; 
this affects the results only at column densities < 10^^ cm~^. We also assume an attenuating 
column of 10^^ cm~^ between the X-ray source and the region being modeled. This value 
is based on the inferred column densities in the cold component (modeled in Section |5]). 
The physical and chemical state of the gas is calculated using an iterative scheme; radiative 
transfer of cooling radiation is handled with an escape probability treatment, including the 
effects of dust trapping of line photons. 

Figure [7] shows XDR models for CO abundance (bottom left) and CO surface brightness 
(bottom right). In the left panel, red lines are the model CO abundances and in the right 
panel red lines are CO surface brightness; in both panels blue lines are the gas temperatures, 
and green lines are the expected/observed values. The black dashed rectangle shows the 
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accepted range of the source size (radius from the center of the XDR) and hmits on cloud 
density (Table 3). It extends up to about 240 pc as calculated from the dust SED in Section 
lU The blue circle is the 1 — a temperature range from the CO non-LTE modeling (Table 
3). Within this temperature range, the XDR models predict a CO abundance of 10~^ - an 
extremely low value. At these X-ray fluxes and gas densities, the high ionization fractions 
result in rapid chemical destruction of the CO, and thus a low equilibrium abundance. In 
addition, the XDR models cannot match the observed CO surface brightness, as shown in 
the bottom right panel. In both cases, much higher densities > 10^'^ cm~^ will be required 
to match the observed values. In fact, if we lower the X-ray luminosity to bring down the 
destruction rate of CO, the temperatures would go down, which is not consistent with the 
very warm temperatures that we derive from the CO emission. So under no condition can 
we produce the observed warm CO at such low densities. This is not evidence against an 
XDR in Arp 220, but shows that an XDR cannot be the source of the warm CO emission. 
Such low density molecular gas can still survive in the presence of an X-ray luminous AGN 
as long as it is shielded from it (from large dust column densities). 

We compare this case with the Cloverlea f quasar, which has a luminous XDR. A similar 
analysis was done by [Bradford et al.l ( 120091 ). in which Z-spec measurements of high- J CO 
lines (up to J = 9-8) combined with low-J lines from other instruments were used to model 
the CO. The density of the gas in Cloverleaf as derived from non-LTE modeling was much 
higher, n ~ 10^ cm~^ and the temperature was much lower, T ~ 100 K. The XDR models 
readily produce the observed CO surface brightness and abundance at the observed density, 
temperature and size scale of the emission. Thus, XDR excitation of CO is possible in 
galaxies, but we have ruled it out as a dominant heating mechanism for the warm molecular 
gas in Arp 220. 



5.3.3. Cosmic Rays 



It is widely believed that the low-energy cosmic ray (CR) protons are efficient in ionizing 
and exciting the gas in the ISM. There is evidence of CR heating of molecular gas in our 
Galaxy (Goldsmith & Langer 1978), and in external galaxies particularly the starbursts such 
as M82 (Suchkov et al. 1993) and NGC 253 (Bradford et al. 2003). 



Meijerink et al.l (120111 ) [MEll] modeled the effect of CR heating on the chemistry of 



molecular gas in the ISM of Arp 220. They give predictions for column densities of CO and 
other species for CR rate s ranging from the G alactic value of 5 x 10~^^ s~^ to 5 x 10~^^ s~^. 



as expected in ULIRGs (jPapadopoulosI |2010[ ) . This range of CR rates are reasonable for 
a ULIRG when compared to the high CR rates of ~ 10^ times the Galactic value, derived 
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from the observations of gamma rays in M82 and NGC 253 flAcero et al.ll2009t lAbdo et al. 



20101). MEll find that the total CO column density is insensitive to CR flux and therefore 



we cannot say much about the effect of CR on the warm CO from these models. However, 
the CR heating rate using the upper limit of the CR rate (5 x 10~^^s~^) from MEll is 
~ 7.5 X 10~^^ergs~^ per molecule, about two orders of magnitude lower than the energy 
budget from H2 cooling of 20 Lq/Mq (~ 1.7 x 10~^^ergs~^ per molecule). In other words, 
a CR rate of ~ 10~^^ s^^ is required to satisfy the total cooling. This r ate is about 
two orders of magnitude above the prediction for ULIRGs (|Papadopoulosll2010l ) . One would 
expect a distribution of CR rates such that the tail of this distribution can provide the higher 
rates required to satisfy the energy budget. However, this would imply that we should also 
see a distribution in temperature of the warm CO gas, which is not reflected in the narrow 
temperature likelihood we obtained from the non-LTE modeling. In summary, CR are not 
likely to have sufficient energy required to explain the observed warm molecular gas. 



5.3.4-- Non-Ionizing Energy Source 



None of the above sources can produce the observed warm CO. We believe that a 
non-ionizing source like mechanical energy is required to heat the molecular gas. Large 
mechanical energy is expected in a merging system like Arp 220. This is also supported 
by a velocity gradient (frms = A^H2^^/^H2) large as 20 km s~^ pc~^ we find from our 
RADEX likelihood analysis. Since we do not know the local size scale corresponding to this 
velocity gradient the amount of mechanical energy cannot be determined. But we can obtain 
a rough estimate of energy from turbulent heating for a region containing 1 of gas of 
density ~ 10^ cm~'^. The length scale corresponding to this volume is about ~ 0.25 pc. Over 
this length scale the turbulent velocity will be about 5 km s~^ (using the velocity gradient 
of 20 km s~^ pc~^). Using this velocity in the expression of turbulent heating per unit mass 
from Bradford et al. J2005h we obtained 0.01 Lq/Mq, which is much lower than the observed 
H2 cooling. Unlike M82 (PIO), in which turbulent heating was likely the major source of 
mechanical energy, it is not the case in Arp 220. 

Finally, we consider the energy fro m supernovae and stel lar winds. Using the observed 
supernova rate of z/sn = 4 ± 2 per year (Lonsdale et al.l 120061 ) and (-Esn = 10^^ erg) in the 
following expression from iMaloneyl (119991 ) 



LsN~3xl0^3(-^)( 



SN 



lyr ^ 10^-^ ergs 



erg s 



(5) 



we obtain total mechanical energy of ~ 100 Lq/Mq (using ~ 3 x 10^ Mq) from super- 
novae in Arp 220. The energy output from stellar winds is similar to supernovae (McCray 
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1987) assuming a Salpeter initial mass function, giving a total of ~ 200 Lq/Mq. This is large 
enough such that a small portion of it will be sufficient to match the observed H2 cooling. 
In conclusion, mechanical energy from stellar winds and supernova are strong candidates for 
heating of the warm molecular gas in Arp 220. 



6. HCN in Absorption: Infrared Pumping 



HCN can be a very good tracer of dense molecular gas (iGao fc SolomonI l2004f ): its 
rotational transitions have critical densities 100 - 1000 times higher than CO. We detect very 
high-J HCN lines, from J = 12-11 to J = 17-16, in absorption, which have critical densities 
around 10^ — 10^° cm~^. Their observed equivalent widths are listed in Table 1, with the 
exception of the J = 13-12 line, which is blended with strong CO and water transitions. Such 
high-J transitions of HCN have never been observed before in other galaxies, in emission or 
absorption. Their low-J counterparts from J = 1-0 to J = 4-3 have been observed in emission 
from ground-based observatories. The left panel of Figure [S] shows the observed HCN ladder. 
The transition from emission to absorption happens somewhere between J = 4-3 and J = 
12-11. 

All the HCN absorption lines fall on the linear part of the curve of growth as shown in 
Figure [HI and therefore their column densit ies can be me asured directly from the observed 



equivalent width per the following relation ( 1Spitzerlll968[ ) 



^ = 8.85 X 10-^3 Ar,A/,fc (6) 

where W\ is the equivalent width, Nj is the column density in the jth level, and / is the 
oscillator strength. The unit of the constant in the above equation is cm. The summed 
observed column density of HCN from the five high-J absorption lines is ~ 5.4 x 10^^ cm~^. 

For estimating the total column density and the physical conditions of the gas that 
this HCN is tracing, we again do a RADEX likelihood analysis, but this time we use the 
column densities of the absorption lines instead of line fiuxes. Using column densities directly 
eliminates any dependence on the source size or the line velocity width. Figure [9] shows the 
results of our likelihood analysis. The temperature and density are not well constrained, so 
we only provide lower limits for them in Table 3. However, the column density is extremely 
well constrained with a total column of A^hcn ~ 2 x 10^^ cm~^. The RADEX analysis 
suggests that the gas is warm, very dense and has a pressure of > 4 x 10^ K cm~'^. 



However, there are several lines of evidence that strongly argue that collisional excitation 
(as assumed in the RADEX models) is not the dominant mechanism for populating these 
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levels: (a) The V2 = 1 level lies 1025 K above the ground. In order for it to be effectively 
populated by collisions, the gas kinetic temperature must be a substantial fraction of this 
value. In addition, the densities must be high, since the vibrational Einstein A-coefficients 
are large (~ 1 km s~^). This will lead to thermalization of the rotational levels at some 
T comparable to 1025 K. However, the rotational temperature of the absorption lines is 
only ~ 300 K. (b) Transitions within the V2 = 1 vibrational state were directly observed by 



Salter et al.l (120081 ) using the Arecibo radio telescope. These also have a low rotational tem- 
perature that is very similar to the ground vibrational state J-levels seen in absorption. This 
again argues that the kinetic temperature cannot be large enough to collisionally populate 
V2=l, and (c) From the observed lines within V2 = 1 levels we can estimate the (v2 = 1, 
J=6) column density. This is 3 times larger than the (v2=0, J=6) column predicted by the 
RADEX collisional excitation analysis, and is a very large fraction of the total V2=0 HCN 
column derived from the RADEX results. This is unphysical. Therefore, we believe that the 
high- J lines are popula ted not by collisional excit ation but instead by radiative pumping of 



infrared (IR) photons (jCarroU fc Goldsmith! Il98ll ). 



The IR pumping mechanism works as follows: Assuming that the collisional excitation 
is negligible, an HCN molecule in v=0, J =0 ground state can absorb a 14 /im photon 
from the dust continuum and be excited to V2=l, J=l state (as per the selection rule of 
A J = ±1). This HCN molecule can decay back either to v=0, J =0 or v=0, J = 2 with 
equal probabilities. If the IR pumping rate is large enough, then the molecule in J = 2 
ground state can absorb another 14 fim photon before decaying to J = 1 and J = 0, and 
can be excited to V2=l, J = 3, which can decay into J = 4 ground state, thereby raising 
the excitation of the ground state rotational levels. This mechanism can be efficient in 
populating the low-J HCN lines. But for an intense radiation field, very high-J ground state 
transitions can also be populated. For IR pumping of the vibrational states to significantly 
affect the population of the ground state rotational J+1 level, the pumping rate out of the 
ground state must exceed the spon taneous radiative decay rate out of the J-l-1 level. This 



gives the following criteria (Eq. 5 of ICarroU Sz GoldsmithI ( 119811 )) for the temperature of the 
radiation field 

/ . Aj+i,j 



eWkTd _ 1 ^ A,^,_^' ^ ' 

where / is the dilution factor to account for both the geometric dilution and optical depth 
effects, Td is the temperature of the radiation field, and A's are the Einstein coefficient for 
spontaneous decay. The hv/k is 1025 K for the V2=l, and substantially larger for the other 
vibrational levels, so it is unlikely that any higher levels will contribute to the IR pumping. 
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To populate J = 17-16, the above criteria will require the temperature of the radiation 
field to be > 350K, assuming / = 1 for a blackbody radiation field. For / < 1 this 
temperature must be even higher. The implication here is that the gas observed in absorption 
in HCN needs to see a > 350 K radiation field. However, this does not mean that we will 
find evidence for this thermal component in our dust continuum data. The line of sight to 
this continuum component is likely to be highly obscured. 

There is direct evidence for IR pumping vi a the 14 nm the H CN absorption feature 
in the Arp 220 spectrum obtained with Spitzer ( iLahuis et al.l 120071 ). The column density 
of this 14 fim absorption is a factor of ~ 15 larger (with 30% uncertainty) than the HCN 
column we obtain from the absorption lines. It is difficult to determine if the gas traced 
by the 14 fim bending-mode absorption feature is the same component seen in the pure 
rotation absorption line, especially given the large optical depths in Arp 220. Additional 
evidence for strong IR pumping comes from the observation of the HNC J = 3-2 line, which 
is overluminous in A rp 220; it is brighter than the HCN J = 3-2 line by a factor of ~ 
2 ( lAalto et al.l 120071 ). Vibrationally excited HCN lines have been also detected in other 
luminous infrared galaxies such as NGC 4418 (jSakamoto et al.ll2010l ). There is also evidence 
for IR pumpi ng of higher-J HCN lines in the host galaxy of the z = 3.91 quasar APM 



08279+5255 f lRiechers et al 



20101: 



WeiB et a. 



work as it is argued by iDownes &: EckartI ( 120071 ) that the dust emission in the western 



20071) . This is particularly relevant to this 



nucleus of Arp 220 is strikingly similar to this quasar. If the excitation of the high-J HCN 
lines in Arp 220 is dominated by IR pumping rather than collisions, then the HCN is not 
necessarily tracing dense gas. 

In summary, we propose that the populations of the high-J levels of HCN seen in ab- 
sorption in Arp 220 are produced dominantly by the effects of radiative pumping, along with 
the additional requirement of intense radiation field with temperature of > 350 K, while the 
low-J levels of HCN (J = 1-0 to J = 4-3 observed from the ground) are affected by both 
collisions and radiative pumping, thereby increasing their column density. To model the 
HCN in Arp 220 we will need to include effects from both IR pumping and collisions in our 
RADEX code. This is beyond the scope of this paper. 



7. Atomic Lines: [C I] and [N II] 

As mentioned in Section [3l we detect two [C I] lines and one [N II] line in emission. The 
[C I] lines are very difficult to observe from the ground and so far they have been detected in 
only a handful of galaxies, including M82 (PIO). The ratio of line strengths (in Jy km s~^) of 
[C I] (2-1) to [C I] (1-0) in Arp 220 after correcting for dust extinction is 1.3 ±0.3, consistent 
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with unity within the error bars. This suggests that these hues are optically thick, which is 
also s upported by the ex citation temperature of 26 K derived from their line temperature 
ratio (IStutzki et al.lll997l ). which is too low compared to the kinetic temperature of the cold 
molecular gas (50 K). For comparison, in M82, using FTS measurements of these lines (PIO) 
we found this ratio to be 2.1 ± 0.2. To calculate the [C I] column density in Arp 220 w e use 
the escape probability radiative transfer models from Figure 2 of IStutzki et al.l ( 19971 ). For 
the observed line temperature (in K km s~^) ratio of [C I] (2-1) to [C I] (1-0) of ~ 0.5, and 
the brightness temperature of ~ 13.5 K of [C I] (1-0) line (assuming 500 km s~^ line width), 
we estimate a column density of Nq ~ 10^° cm^^. The low excitation temperature of C I 
suggests that it is associated with the cold molecular gas traced by low-J CO lines with the 
ratio of Nc/Nqo ~ 1. 

The ratio Nq/Nqo is much higher in Arp 220 compared to the Galactic PDRs, e.g, this 
ratio is around 0. 17 for the Or i on bar (Tau ber et al. 1995). Previous studies of atomic carbon 
in galaxies (e.g., IWilsonlll997l : Ilsraelll2005l ) suggest that [C I] emission gets stronger for more 
luminous systems with starbursts and AGN. For com parison, the r a tio of N[c]/Nco in M82 
and the Cloverleaf Quasar are 0.5 and 1, respectively. iKeene et al.l (119961 ) found one-to-one 
correlation between [C I] and ^^CO line strengths for Galactic PDRs. They argue that simple 
PDR models in which the cloud has a plane parallel geometry illuminated on one side by UV 
field cannot produce this correlation. In these simplified models, [C I] forms a narrow layer 
between [C II] on the outside and CO deeper in the molecular cloud, with column density 
of about 5 X 10^^ cm~^. The column density and line strength of [C I] (1-0) line in these 
models are relatively insensitive to physical parameters such as temperature, density and 
the radiation field. The observations of some Galactic PDRs have shown that the [C I] can 
extend f ar deeper in to the molecular clouds, much further than predicted by simple PDR 
models. Ilsraell ( 120051 ) found from the survey of atomic carbon in nearby luminous galaxies, 
that the correlation between [C I] and ^^CO still holds but it is not one to one anymore. 
Instead [C I] emission gets stronger compared to CO for more luminous systems. This is 
shown in their Figure 2 in which the ratio of [C I] (1-0) and ^'^CO (2-1) is plotted as a 
function of [C I] luminosity. On this plot, Arp 220 falls on the upper right corner (with 
coordinates ~ (3.2,4.5)), consistent with the correlation. 

The extinction corrected luminosity of the [N II] 205 fim line is ~ 2.96 x 10^ Lq, which is 
10 times higher than that of M82 (PIO). But compared to their respective FIR luminosities, 
the L[Nii]/i^FiR ratio of ~ 1.7 x 10^^ in Arp 220 is lower than M82 by a factor of 3. The deficit 
of [N II] 205 fim line relative to LpiR in Arp 220 is surprising because [N II] is produced in 
the H II regions and not in the photo-dissociation re gions, where the [C II] 158 fim deficiency 
is generally found in ULIRGs ( iMalhotra et al.lll997l ). However, the [N II] line deficit we find 
here is consistent with the recent FIR line deficit results from Herschel- PACS presented in 
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Gracia-Carpio et al.l (120 111 ). These authors found that the [N II] 122 fim line shows deficits 



similar to the [C II] 158 fim line deficit relative to LpiR for a large sample of galaxies. They 
claim that this is the case with all FIR lines but their case for an [OIII] line deficit is weak 
compared to [N II] line deficit. Their modeling results suggest that the line deficit in H II 
regions is the result of a higher ionization parameter in ULIRGs compared to normal or 
starburst galaxies. A higher ionization parameter will increase the dust absorption of the 
UV photons, which will decrease the fraction of photons available for ionizing the neutral 
gas thereby decre asing the fluxes of [N II] and other FIR lines relative to the continuum 
kbel et allboooh . 

The line ratio of [N II] 122 fim to 205 fim is an excellent density probe of the ionized 
gas. However, in A rp 220 the [N II] 122 /im appears in absorption in the ISO spectrum 
( iBrauher et al.ll2008[ ). so we cannot use this line to estimate the electron density. Instead we 
com pare the line strength of the [N II] 205 line to the [C II] 158 /im line measured by 
ISO (IBrauher et al.ll2008l ). The line ratio of [C II] /[N II] gives an estirnate of the fraction of 
[C II] in the ionized medium compared to the PDRs (e.g., lOberst et al.l ( 120061 )). We find this 
line ratio (co r rected for dust extinction) to be 17 ± 3. We compare this ratio to Figure 2 of 



Oberst et al.l (120061 ). in which they plot the model line intensity ratio as a function of ionized 



gas density. Because we cannot estimate the ionized gas density in Arp 220, we consider the 
minimum and maximum values of this line ratio over the entire density range. These model 
line ratios range from 2.5 - 4.3. Comparing this model line ratio to our observed value of 
17, we constrain the contribution from H II regions to the observed [C II] 158 fim emission 
to be between 15% and 26%. This sugge sts that the naajority of [C II] arises in the PDRs, 
as concluded in the previous works (e.g., iMadden et al.lll993l : ISauty et al.lll998l ). 



8. The Molecular Outflow and Water Chemistry 
8.1. Molecular Outflow: Observations of P Cygni Profiles 



In the standard scenario of galaxy evolution, some of the fundamental open questions 
are: How do feedback mechanisms regulate star formation and supermassive black hole 
growth to pr oduce the observed M — a relation betweeii stell ar velocity dispersions and black 
hole masses (IFerrarese fc Merrittll2000l : iGebhardt et al.ll2000l )? Can the energy injection from 
starbursts or AGN quench star formation by expelling the gas from the host galaxy? Is this 
the process by which massive gas-rich galaxies deplete their interstellar media to become 
gas-poor early-type galaxies? Direct observational evidence of massive molecular outflows 
associated with AGN/starbursts can significantly impact our understanding of the connection 
between feedback and galaxy evolution. 
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Observations of nearby galaxies with AGN have shown outflows but generally only in the 
ionized gas component, which is a minor fraction of the gas mass in the host galaxy. Direct 
observational evidence of massive molecular outflows was recently discoyered in two galaxies 



using the IRAM PdBI and Hersche l PAC S: Mrk 231 flFischer et al 



20 



a 



Feruglio et al. 



2010h detected 



201(1 ) and NGC 1266 flAlatalo et al.ll201lh . In Mrk 231, iFischer et all | 
massive outflow with velocity of ~ 1400 km s~^ using the P Cygni profile seen in the OH 



lines. Independently, iFeruglio et al.l (120101 ) also detected the same outflow using the CO 
(1-0) line, and derived an outflow rate of 700 Mfn /year, which will dep lete the gas reservoir 
in 10^ yrs. In NGC 1266, a nearby field SO galaxy, lAlatalo et al.l (l201l[ ) discovered a massive 
molecular outflow originating from the nucleus. They found the outflow velocity of 400 km 
s~^ and outflow rate of 13 Mq/jt with a depletion time scale of 85 Myr. They conclude that 
the outflow has to be drive by an AGN in absence of sufficient star-formation. 

A P Cygni line profile (redshifted emission peak with corresponding blue-shifted absorp- 
tion) is a strong signature of outflow activity. We detect P Cygni profiles in the OH"*", 1120"*" 
and HF lines as shown in Figure [TOl suggesting molecular outflow. Since our P Cygni line 
profile is not resolved we can only put an upper limit on the outflow velocity of about 250 km 
There is independen t evidence for molecular outflow from ground-based observations 



by ISakamoto et al.l (l2009l ). in which P Cygni profiles were detected in HCO"'" lines. They 
found that the outflow has a velocity of 100 — 200 km s^^ and is confined within a radius of 
50 pc. The total HCO^ column density was about 10^^ cm~^. The outflow in HCO^ is very 
likely the same outflow that we detect in OH"*" and H20"^. Our upper limit of 250 km s~^ 
on the outflow velocity is a limit on the line-of-sight velocity. It is likely that the outflow 
has a wide angle and is not in the form of a narrow jet. In this situation, the maximum 
outflow velocity will be higher than 2 50 km s~^ dep e nding on the inclination of the disk. 
Using the inclination angle of 45° from IScoville et al.l (119971 ) and assuming that the outflow 
is orthogonal to the disk, the maximum outflow velocity is ~ 350 km s~^. As this is lower 
than the escape velocity of the nuclear region of Arp 220 (~ 420 km s~^ at a nuclear radius 
of 250 pc), we believe the outflow is gravitationally bound to the nuclear region. 



8.2. Importance of H20^ Detection 

The detection of H2O"*" by Herschel in the Galactic and extragalactic ISM was very 
surprising. This is because H2O"*" is highly reactive and is expected to rapidly react with H2 
to form HaO"*". The absence of HsO"*" and a presence of strong absorption of H2O"*" requires 
a gas with a low molecular fraction, which can be the result of high cosmic ray or X-ray 
fluxes. The detection of H20"^ and its implication on the physical conditions of the molecular 
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gas is reviewed by Ivan Dishoeck et al.l (1201 if ) . The presence of strong absorption features in 
both 9^20"*" and OH"*" also suggest that the production of water in the gas associated with 
the outflow is dominated by low temperature ion-neutral chemistry. At high temperatures 
(T > 500 K), the water production is dominated by neutral-neutral chemical reaction (no 
production of H2O+ and 0H+), which increases the abundance of water significantly. 



8.3. Origin of the Ionic Molecular Species: Evidence for AGN in Arp 220 



The presence of a hidden AGN in Arp 220 has long b een debated. The X-r ay obser- 
vations from ROSAT (IHeckman et al.l Il996l ) and Chandra (jClements et al.l |2002|) reveal a 



relatively weak X-ray source with a luminosity of 4 x 10 ergs s . Neither study can rule 
out a hidden AGN that is obscured by the la rge column of N('H2) 



10^^ cm ^ found in 



this study and in previous observations (e.g., iGonzalez- Alfonso et al.ll2004j ). In addition, 
Haas et al.l ( 1200 ll ) found an exceptionally low ratio of 7.7 /xm PAH flux to submm contin- 
uum for Ar p 220 and argue that a hi dden AGN is required to power much of the luminosity 
of Arp 220. iDownes &: EckartI (120071 ) found evidence for a hot, compact, and optically thick 
dust ring in the western nucleus from 1.3 mm I RAM data, which t hey interpreted as being 
heated by an AGN accretion disk. More recently, lEngel et al.l ( 120111 ) found that this putative 
AGN's position is coincident with the stellar and CO (2-1) kinematic center, suggesting that 
the stars in the western nucleus are gravitationally bound to either a supermassive black 
hole or an extremely dense, young nuclear starburst. 

Here we consider both CR and XDR models to determine which one of these sources can 
match the observed column densities of the ions involved in the chemical network producing 
water. The absorption lines of these ionic species are saturated; their equivalent widths are 
not linear on the curve of growth. They appear saturated because of the lower resolution of 
our FTS spectrum. Hence, we can only determine the lower limits on their column densities 
from their equivalent widths (listed in Table 2) using the relation in Equation [61 We obtain 



Ar(OH+) > 10 



15.6 2 



A^(H20+) > lO^^-i cm-2, and N{R20) > 14.3 cm'^. These lower 



limits are calculated from the ground state transitions of these molecules, which contain 
a majority of the population. The HaO"*" line is absent in our spectrum so we calculate 
a 3—0" upper limit of N(H30"'') < 10^^ '^ cm~^. We co r apare d the lower limits, A^(OH+) 
and iV(H20"'"), with the CR models of iMeijerink et al.l (120111 ) (see their Figure 7). They 



explored a wide range of the CR ionization rates as discussed in 15.3.31 We find that none of 
their models (with low and high densities; low and high Go fields) can produce the observed 
column densities for any value of CR ionization rates. These CR models underestimate these 
column densities of OH"^ and H20^ by almost an order of magnitude, and overestimate the 
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H3O"'" column. 

We now consider XDR models, which were explained in 15.3.21 In Figures [TT] a - e, we 
present predictions for the column density over the nn and A^h grid for all the ions including 
HCO"'" that are involved in the outflow. The allowed regions of columns densities are shown 
by shaded area. This model was run for a source size of 50 pc (estimated for HCO"'" by 



Sakamoto et al.l (120091 )). In Figure [TTf. total column densities of each species is shown as 
a function of total hydrogen column and the observed lower limits are marked by dashed 
line^. Compared to CR, the XDR models can readily produce the column densities of all 
these species. In order to produce the observed columns, the luminosity of the XDR has to 
be ~ 10"^^ ergs s~^. 

In figures [TT] a - e, we note that OH"'' column density cannot be produced for rin > 10^ 
cm~^. This gives us an upper limit on the density of the gas associated with the outflow. 
This limit will change with source size. We tested this by running XDR models for 100 pc 
and 240 pc source sizes. With a larger source size the density limit after which OH"^ column 
density cannot be matched goes down. It is also possible that there is a density gradient in 
the outflow, in which case the HCO"*" emission could be dominated by a denser component. 
There is support for this argument in Figure fTTT . where the column densities of all the ions 
except HCO"*" correspond to the hydrogen column density of ~ 2 x 10^^ cm~^. The HCO"*" 
column density corresponds to a much higher column density of ~ 6 x 10^^ cm~^ suggesting 
that HCO"^ could be coming from denser part (nn ~ 10^ cm~'^) of the outflow. 

We use the hydrogen column density of ~ 2 x 10^^ at which the column densities of 
all the ions (except HCO"*") are satisfied, to estimate the lower limit on the amount of mass 
associated with the bound outflow. We find this lower limit to be Moutflow ^ 2 x 10^ M©, 
using a source size of 50 pc and conservatively assuming that the outflow only covers the 
continuum source. 

One of the interesting things about the Arp 220 outflow is that it appears to be very 
different from Mrk 231 and NGC 1266. This is because the mass associated with the outflow 
in Arp 220 is large but it will not escape the nuclear region of this galaxy. Does this mean 
that the standard scenario of AGN/starbursts driving the molecular gas out of the galaxy 
does not apply here or is Arp 220 in a different evolutionary stage compared to Mrk 231 
such that at some time in the future this gas will be expelled out of the galaxy? 

In conclusion, the CR models not only fail to produce the observed column densities of 
H2O"'" and OH"^ by almost an order of magnitude but they also over estimate HsO"*" by a 



^For H3O"'" and HCO"'", the dashed hne represents an upper hmit and a total column density respectively. 
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large amount. An XDR with luminosity of 10^^ ergs s ^ can easily match the observations 
of all the above ionic species, strongly favoring the existence of an AGN in Arp 220. 



9. Discussion: Is Arp 220 a Good Template for High-z Submillimeter 

Galaxies? 



Arp 220's spectrum is commonly used as a template for high-z submm galaxies be- 
cause it is bright, and therefore has been extensively studied across the electromagnetic 
spectrum. We discuss here in the light of new results from this study on dust and molec- 
ular gas, whether Arp 220 is still a good template. We compare some characteristics of 
Arp 220 with several high-z SMGs for which the molecular and co ntinuum data were re - 
cently published. There are four high-z sour ces, SMMJ2135-0102 (jPanielson et al.ll201ll ). 
GN20 jcarilh et al.l boioh. SMM18423+5938 Jlestrade et aPboioh . and SMM10571+5730 
( IScott et al.ll201ll ). with z ranging from 2.3 - 4.0, for which submm spectroscopic data was 
obtained using several ground-based facilities (Z-Spec, IRAM-PdBI, APEX, SMA, CARMA, 
etc.). Several CO transitions from low- J to mid- J were observed for all four systems. In all 
cases, the CO SLED peaked around CO (6-5) and then turned over, similar to Arp 220 
and M82. Does this mean that the CO excitation in these high-z systems is dominated by 
starbursts an d not AGN? For the latter, the CO SLED is not expected to turn over aroun d 



mid- J levels ( Ivan der Werf fc Spaand l2009t [Bradford et al.l 120091 : Ivan der Werf et al 



More recently, EV LA observations of five high-z quasar host galaxies (IRiechers et al. 



20101). 



2011 



Ivison et al.l 1201 if ) indicate that SMGs with AGN show no cold gas component; i.e., the CO 
J = 1-0 line strengths in these five systems are consistent with warm, highly excited gas 
that is also seen in the higher-J CO lines. They conclude that the gas-rich quasars and 
AGN-free SMGs represent different stages in the early evolution of massive galaxies. It is 
complicated to compare this scenario directly with our findings in Arp 220, in which the 
CO SLED shape is consistent with gas that is excited by star-f ormation bu t obser vations of 
other molecules also indicate the presence of a luminous AGN. iLupu et al.l ( 120101 ) observed 
five submm-bright lensed sources recently detected by the Herschel Astrophysical Terahertz 
Large Area Survey (H-ATLAS) using Z-Spec. Up to 3 mid-J CO transitions were detected 
in these systems, not enough to compare the shape of the SLED to Arp 220. 

We compare the ratio of -Z^co,j=6-5/-^fir in Arp 220 with a sample (discussed above) 
of 9 high-z sources (using individual CO h ne luminosities) and find that this r atio does not 
vary much: it is ~ a few xlO~^. In six ( IScott et al.l 1201 ll : iLupu et al.l |2010[ ) out of nine 
galaxies, the modeling of CO was similar to this work, but the results from this modeling 
have large error bars because fewer lines were available. This makes it hard to conclude 
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if all these higl i-^: systems also h ave warm molecular gas as observed in Arp 220. The 
CO analysis by IScott et al.l (120 111 ), for which six CO transitions were available, indicates 
that HLSW-01 has significant amount of warm molecular gas; the CO SLED was fit well 
by single temperature component with Tkin ranging from 86 - 230 K, with no evidence for 
a cold temperature component. The data in this system are consistent with no cold gas 
component. In addition, we compare th e ratio of Lct/-^ftr ~ 6.7 x 10~^ in Arp 220 with 
a recent survey of CI in high- 2; galaxies ( Walter et al.l |20 111 ). This study found the ratio to 
be ~ 7.7 X 10^^ from a sample of high-redshift galaxies. The two ratios are consistent with 
each other and are roughly equal to the luminosity of a single CO line in Arp 220, implying 
that CI is not an important coolant compared to CO in Arp 220 and the high-z sample. 



In the six (jScott et al.l l2011t iLupu et al.l I2OIOI ) out of nine galaxies, the modeling of 
dust was similar to this work; i.e., they were modeled using Equation 1. All these six high-^ 
systems show warm dust emission with Tdust varying from 55 K - 90 K, and high optical 
depths of Td ~ 3.5 — 5 at 100 /xm. However, for most of these galaxies (3 was fixed at 2.0; 
since (3 is correlated with Tdust (as shown in Figure H]), the dust temperatures could be lower. 
Both the dust temperature and high optical depth in Arp 220 are consistent with the range 
of temperature and optical depth estimated for the high- 2; systems. Due to large optical 
depths, dust extinction correction should be considered for the CO lines observed in these 
SMGs. As demonstrated in this work, the kinetic temperature of the warm molecular gas is 
grossly underestimated (see section 5.2) if the CO line fluxes are not corrected for extinction. 
We note that this high-z sample is small and suffers from a selection effect which favors the 
detection of warm galaxies for a given dust mass. 

In summary, even though FTS observations of Arp 220 indicate that it is a peculiar 
ULIRG with unusually large dust optical depths, massive molecular outflow, strong IR pump- 
ing and warm molecular gas, the global properties such as CO line to FIR continuum ratios, 
dust temperatures and shape of the CO SLED seems to be consistent with a sample of high- 2; 
SMGs. Currently, the SMG sample is small and the S/N of the line observations is not as 
good as for Arp 220. It will be important to revisit this comparison as more molecular line 
observations of dusty, luminous, high-z become available. 



10. Conclusions 

In this paper, we presented a complete spectral view of Arp 220 from 190 - 670 /im using 
the SPIRE-FTS; most of this wavelength region has not been accessible prior to Herschel. 
The continuous wavelength coverage provided by the FTS eliminates any systematic uncer- 
tainties due to cross-calibration issues arising from using multiple instruments observing over 
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different wavelength windows, as has been the case for ground-based observations of nearby 
galaxies for a long time. The wavelength region covered by the FTS contains the higher-J 
rotational levels of CO, HCN and other molecules, which were crucial for carrying out a 
comprehensive characterization of the molecular gas. The continuous wavelength coverage 
also facilitated detections of rare molecular and atomic species that a targeted narrow band 
search would not have achieved. In addition to CO and water species, we report unexpected 
detections of very high- J absorption lines of HCN. The spectrum also show strong absorption 
in 0H+, H2O''", CH+, HF, and several nitrogen hydrides. The two emission features of H2O+ 
at 742 and 746 GHz have been seen for the first time. We also detected C I and [N II] lines. 
These atomic and molecular species are very difficult to observe from the ground, and in 
some cases completely inaccessible. In addition to line detections, the FTS also provided 
a good measurement of the continuum, which is invaluable in characterizing dust in the 
nearby starbursts and ULlRGs. Hence, FTS instrument is playing a very important role 
in providing information at wavelengths that will be complementary to ALMA and other 
ground- and space-based telescopes. 

Below we list the conclusions for Arp 220, which are drawn from the analysis of the 
FTS spectrum: 

(a) Modehng the continuum revealed dust that is warm with T = 66 K, and has an un- 
usually high optical depth, with Td^ioofim. ~ 5. We derive a dust mass of about 10^ Mq and a 
high total hydrogen column density of ~ 10^^ cm~^ with an uncertainty of a factor of ~ 3, 
completely dominated by the uncertainty in the dust cross section. 

(b) Warm Molecular gas: The extinction corrected observed CO luminosity is dominated 
by the mid-J to high-J lines, peaking around CO(6-5). The non-LTE radiative transfer 
modeling strongly indicates that the mid-J to high-J lines are tracing warm molecular gas 
with Tkin~ 1350 K. The low-J transitions arc tracing cold gas with Tkin~ 50 K. The inferred 
temperature for the warm and cold components are much lower if CO line fluxes are not 
corrected for dust extinction. These two components are not in pressure equilibrium and 
the observations show that their line widths are different by a factor of ~ 1.5. The mass 
of the warm gas is about 10% of the cold molecular gas but dominates the luminosity as 
well as the cooling over the cold CO. The ratio Lco/-^fir ~ 10~^, where Leo is the total 
CO luminosity. The temperature of the warm molecular gas is in excellent agreement with 
the temperature derived from H2 rotational lines observed from Spitzer, implying that CO 
is still a good tracer of H2 at these high temperatures. At 1350 K, H2 dominates the cooling 
of the ISM over CO. We also conclude that the contribution of the dense gas to the observed 
CO emission is small. 
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We have ruled out PDRs, XDRs and cosmic rays as possible sources of this warm molec- 
ular gas within the context of models we compared our data to. The mechanical energy from 
supernova and stellar winds can satisfy the energy budget required to heat this gas but we 
still do not know the exact mechanism that hea ts this gas. Such wa rm molecular gas has 
been confirmed in only two galaxies so far, M82 (jPanuzzo et al.ll2010l ) and Arp 220, as part 
of the submm wavelength region covering high-J transitions of CO was made accessible by 
SPIRE. In both cases some type of non-ionizing source is required to heat this gas. As FTS 
data for more galaxies become available, this issue will be investigated further. 



(c) The very high-J lines of HCN appear in absorption. The transitions from emission 
to absorption takes places somewhere between J = 4-5 and J = 12-11. These high-J lines are 
populated via IR pumping of photons at 14 /xm. The condition for IR pumping to populate 
J = 17-16 level requires an intense radiation field with T > 350 K. 

(d) Massive molecular outflow: The signature of a molecular outflow is seen in the P Cygni 
profiles of 0H+, H2O"'" and H2O; major molecules involved in the ion-neutral chemistry pro- 
ducing water in the ISM. The outflow has a mass of ~ lO'^ M0 and velocity < 250 km s"^ 
It is massive but bound because its velocity is less than the escape velocity of the Arp 220 
nuclei. The upcoming Atacama Large Millimeter Array (ALMA) will have the capability to 
map this outflow. A high resolution follow-up from HIFI and ALMA is required to accu- 
rately characterize this outflow. So far, massive molecular outflows have only been detected 
in Mrk 231 and NGC 1266. 



(e) The existence of an AGN in Arp 220 has long been debated. We found significant 
evidence for an AGN in Arp 220. Our modeling shows that the large observed column den- 
sities in OH"*", H2O"'' and H2O can only be produced by a luminous XDR with Lx = 10^^ 
ergs s^^. The outflow can be associated with either this AGN or the starburst. 

(f) The [N II] 205 /im line shows a deficit relative to the FIR continuum, similar to the 
[C II] deficit found in ULIRGs. This is surprising because [N II] arises in the H II regions 
and not in PDRs. This deficit is consistent with recent results from a survey conducted by 
Herschel-PACS, in which deficits were found for the [N II] 122 /xm line for several ULIRGS. 
The line deficit in the H II regions could be a result of higher ionization parameter in ULIRGs 
compared to more quiescent systems. The observed line ratio of [C II] /[N II] constrains the 
percentage of [C II] emission in PDRs compared to H II regions to be between 85% - 75%. 

(g) The ratio of Nq/Nqo is close to unity, which is high compared to the Galactic PDRs, but 
consistent with the correlation between [C I] /CO line strength and FIR luminosity found 
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using data from several nearby luminous galaxies with starbursts and AGN. 

In summary, this is the first time we have a comprehensive picture of the state of molecular 
gas in Arp 220. We find that the molecular gas is highly influenced by the mechanical energy 
provided by the on-going merger. It is important to learn more about the geometry of this 
galaxy to see how the warm gas and the outflow are distributed. Observations from ALMA 
will be able to provide this information. 
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Fig. 1. — Herschel SPIRE-FTS spectrum of Arp 220. The spectrum shows the FIR contin- 
uum between 190 - 670 /xm in (a). Red sohd points in (a) are the continuum measurements 
from SPIRE photometer and the dotted curves show the photometer bandpasses with arbi- 
trary normalization. Line identifications are shown for several molecular and atomic species 
in (b) and (c) with like colors for like species. 
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2. — A zoom- in view of some of the spectral lines detected in the Herschel-SPIRE 
rum of Arp 220. 
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Fig. 3. — Dust spectral energy distribution of Arp 220: The dot-dashed hne is a modified 
blackbody model fit to the combined FTS, SPIRE-photometer and ISO-LWS data. Other 
data from the literature are over plotted for comparison and agree well with the fit. A 
single temperature component model fits the data extremely well down to 15 /im, below 
which it breaks down. The size of the data points is approximately equal to the size of the 
uncertainties (random and systematic). 
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Fig. 4. — Modeling the dust SED yields tight constraints on a single- component dust spec- 
trum in Arp 220. The contours show 68%, 95% and 98% confidence levels. 
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Fig. 5. — Left: Extinction corrected luminosity distribution of the CO ladder from J = 1-0 
to J = 13-12 with the exception of J = 10-9 line, which is blended with a water line. The 
black solid circles are FTS measurements and the blue triangles are average line fluxes from 
ground-based measurements. Also shown for comparison are new measurements by Z-spec 
and JCMT. Right: Non-LTE model fit to the observed CO line temperatures. There are two 
temperature components: a warm component (dotted line) traced by the mid-J to high-J 
lines and a cold component (dashed line) traced by the low-J lines. 



•a 
o 
o 

J3 



o 



1.2 

1.0 

0.8 

0.6 

0.4 

0.2 
0.0 



10 






1000 



2.0 2.5 3.0 3.5 4.0 4.5 5.0 18.50 19.25 20.00 20.75 21.50 



Log n[H2] (cm ^) 



LogNco(cm ) 




Log P (K cm"') 



Fig. 6. — Radiative transfer modeling of CO: Likelihood distributions for gas kinetic tem- 
perature, density, column density and pressure. The blue line represents the cold component 
at 50 K obtained from modeling low-J transitions and the red line represents the warm 
component at 1350 K obtained from modeling the mid-J to high-J transitions. 
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Fig. 7. — Top: PDR models for CO ( iWolfire et al.ll2010h . The red lines are model contours, 
green line is the observed value (except for the CO abundance for which it represents an 
expected value derived from Galactic observations), and the black dashed lines are acceptable 
density limits from non-LTE modeling of CO. Bottom: XDR models for CO abundance and 
CO surface Brightness in erg cm~^ s~^ sr~^. The red lines are model contours, blue lines are 
temperature contours, green lines are observed values, black dashed region shows acceptable 
ranges for the source size and density from Table 3, and the blue dash circle shows the 3 — cr 
limits for the temperature as listed in Table 3. 
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Fig. 8. — Left - HCN rotational ladder: The low-J lines measured from the ground are in 
emission, while the high-J lines measured by SPIRE-FTS are in absorption. Right - HCN 
curve of growth: All the HCN absorption lines from J = 12-11 to J = 17-16 are on the linear 
part of the curve of growth. 
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Fig. 9. — Radiative transfer modeling of HCN: Likelihoods for gas kinetic temperature, 
density, column density and pressure. These results do not include the effects of IR pumping. 
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Fig. 10. — P Cygni detected in 0H+, H20+ and HF lines, suggesting molecular outflow in 
Arp 220. The dashed region highlights the P Cygni profile of the molecule. 
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Fig. 11. — XDR Model: (a) - (e) Predictions for columns densities of several ion species 
over the nn and A^h grid. The shaded area shows the allowed regions of column densities 
for 0H+, H20"^ and H2O HsO"^ and HCO"*". (f) - Total columns densities of each species 
as a function of total hydrogen column. The dashed lines are observed lower limits on the 
column densities for all molecules except for HCO^ and HaO"*", for which they represent an 
observed value and a 3 — a upper limit, respectively. 



-52 - 



Table 1. Line Strengths of Emission and Absorption Lines in Arp 220 



Emission Lines 



ID 


Transition 


Rest Frequency 


Flux 
(Jy km s~l) 


Source 


Comments 


i^CO 


J = 4 - 3 


461.04077 


4550 ± 330 


JPL/CDMS 




12 CO 


J = 5 - 4 


576.26793 


3660 ± 170 


JPL/CDMS 




12 CO 


J = 6 - 5 


691 47,308 


4070 -1- 80 


JPL/CDMS 




12CO 


J = 7 - 6 


806.65181 


3460 ±170 


JPL/CDMS 




12CO 


J = 8 - 7 


921.79970 


3260 ± 200 


JPL/CDMS 




12CO 


J = 9 - 8 


1036 91239 


2920 -1- 270 


JPL/CDMS 




12 CO 


J = 10 - 9 


1151.98545 




JPL/CDMS 


blended with H2O 3i2 — 221 


12 CO 


J = 11-10 


1267.01449 


1280 ± 90 


JPL/CDMS 




12 CO 


J = 12-11 


1381.99510 


840 ± 60 


JPL/CDMS 




12 CO 


J = 13-12 


1496.92291 


450 ± 100 


JPL/CDMS 




CI 


J = 1-0 


492.16065 


1840 ± 70 


CDMS 




CI 


J = 2-1 


809.34197 


2130 ± 400 


CDMS 




H20 


2ii — 2o2 


752.0332 


2970 ± 110 


LOVAS 




H20 


2o2 — 111 


987.92676 


3440 ± 450 


JPL 




H20 


3l2 — 3o3 


1097.36479 


3080 ± 160 


JPL 




H20 


3l2 — 221 


1153.12682 


3460 ± 330 


JPL 


blended with CO J = 10-9 line 


H20 


321 — 3l2 


1162.91159 


4930 ± 300 


JPL 




H20 


422 — 4l3 


1207.63871 


1750 ± 170 


JPL 




H20 


220 — 2ii 


1228.78877 


1980 ± 300 


JPL 




H20 


523 - 5l4 


1410.61807 


940 ± 110 


JPL 




H20+ 


2o2 - 111, ^5/3-3/2 


742.0332 


840 ± 90 


BRlOb 




H20+ 


2o2 - 111, J3/2-3/2 


746.1938 


980 ± 130 


BRIO/JPL 




HCN 


J = 6-5 


531.71635 


1700 ± 120 


CDMS 


(?) 


HCO+/HOC+ 


J = 6-5 


536.82796 


850 ± 160 


CDMS 


(?) 


Nil 


3Pi -3 Po 


1461.1319 


1800 ± 120 


C93<= 


very broad 


Absorption lines 


ID 


Transition 


Rest Frequency 
(GHz) 


(xl0-2;im) 


Source 


Comments 


HCN 


J = 12-11 


1062.980 


2.5 ± 0.2 


CDMS 




HCN 


J = 13-12 


1151.449 




CDMS 


blended with CO (10-9) and water line 


HCN 


J = 14-13 


1239.890 


1.9 ± 0.2 


CDMS 




HCN 


J = 15-14 


1328.302 


1.6 ± 0.1 


CDMS 




HCN 


J = 16-15 


1416.683 


1.2 ± 0.1 


CDMS 




HCN 


J = 17-16 


1505.030 


0.9 ± 0.1 


CDMS 




CH+ 


1 - 


835.07895 


18.8 ± 0.3 


CDMS 




OH+ 


l()i — O12 


909.1588 


11.1 ± 1.4 


CDMS 




OH+ 


I22 — Oil 


971.8053 


16.6 ± 2.0 


CDMS 




OH+ 


I12 — 0l,2 


1033.118 


16.5 ± 1.8 


CDMS 




H2O 


lio — Ooo 


1113.34296 


8.4 ± 1.4 


JPL 




H2O+ 


111 - Ooo, J3/2-1/2 


1115.2040 


9.1 ± 1.2 


GIO^ 





Table 1 — Continued 



Emission Lines 



ID Transition Rest Frequency Flux Source Comments 

(GHz) (Jykms-l) 



H2O+ 111 - Ooo, ^3=1/2-1/2 1139.5606 8.6 ± 0.5 GIO 

HF J = 1-0 1232.47622 7.0 ± 0.7 JPL 



''JPL: Jet Propulsion Lab, CDMS: Cologne Database for Molecular Spectroscopy, LOVAS: NIST recommended rest frequen- 
cies. 

^Bruderer et"aL I l l2010h 
teupta et ai] ll2010h 

■ iBrown et al] lll994l' ): IColgan et al.l lll993l) 

"The line fluxes can be corrected for dust extinction using the following two relations: = {u/uq)^ (where /3 = 1.84 and 
VQ = 1270 GHz), and the mixed dust extinction model given by / = /o(l — e~^'*)/''"[j. The extinction correction factor {Iq/I) 
ranges from 1.076 at 450 GHz to 1.953 at 1600 GHz 



Table 2. Continuum fluxes 



SPIRE-FTS 



Wavelength 


Flux 


(/im) 


(Jy) 


586 ± 70 


2.2 ± 0.5 


480 ± 46 


4.3 ±0.9 


400 ± 32 


7.6 ± 1.2 


345 ± 24 


12.2 ± 1.6 


300 ± 18 


18.4 ±2.0 


270 ± 14 


24.4 ±2.2 


240 ± 12 


32.2 ± 2.4 


220 ± 10 


39.7 ±2.4 


200 ±8 


50.9 ± 1.6 



SPIRE-photometer 



250 
350 
500 



30.1 ± 1.5 
11.7 ±0.6 
3.9 ±0.2 
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Table 3. Highest Likelihood Values and Their Uncertainties from Non-LTE Modeling 



CO - Warm 



Quantity 


Highest Likelihood Value 


1— (T range 


3— (T range 


Tkin (K) 


1343 


1247 - 1624 


1174 - 1699 


Logio"(H2) (cm 8) 


3.2 


3.0 - 3.2 


3.0 - 3.3 


LogioAf(i2CO) (cm-2) 


19.4 


19.4 - 19.5 


19.4 - 19.6 


LogioP (K cm-3) 


6.3 


6.2 - 6.4 


6.1 - 6.4 


Mgas (xIO^Mq) 


4.7 


4.6 - 6.0 


4.3 - 6.9 


Lco(x10«Lq) 


2.0 


1.9-2.1 


1.7 - 2.3 


dv/dr (km s~^pc~^) 


20 


10 - 24 


8-28 


CO - Cold 




50 


34 - 67 


32 - 83 


Logion(H2) (cm ^) 


2.8 


2.6 - 3.2 


2.6 - 3.40 


LogioJV(i2CO) (cm-2) 


20.3 


19.9 - 20.3 


19.8 - 20.4 


LogioP (K cm^'^) 


4.5 


4.5 - 4.8 


4.4 -- 5.0 


Mgas (xIO^Mq) 


5.2 


2.0 - 5.4 


1.5 - 6.8 


dv/dr (km s~^pc~^) 


1.4 


1.3 - 6.0 


1.1 - 8.5 


HON - high-J 


Tkin (K) 


> 320 






Logion(H2) (cm ^) 


S 6.3 






LogioiV(i2CO) (cm-2) 


15.3 


15.3 - 15.4 


15.2 - 15.5 


LogioP (K cm-3) 


>9-6 







